Searches for High Frequency Variations 
in the 8 B Solar Neutrino Flux at the Sudbury Neutrino 

Observatory 

B. Aharmim 6 , S.N. Ahmed 14 , A.E. Anthony 17 -", N. Barros 8 , E.W. Beier 13 , A. Bellerive 4 , 
B. Beltran 1 , M. Bergevin 7 - 5 , S.D. Biller 12 , K. Boudjemline 4 , M.G. Boulay 14 , T.H. Burritt 19 , 
B. Cai 14 , Y.D. Chan 7 , D. Chauhan 6 , M. Chen 14 , B.T. Cleveland 12 , G.A. Cox 19 , 
X. Dai 14 - 12 - 4 , H. Deng 13 , J. Detwiler 7 , M. DiMarco 14 , P.J. Doe 19 , G. Doucas 12 , 
P.-L. Drouin 4 , C.A. Duba 19 , F.A. Duncan 16 - 14 , M. Dunford 13 - 6 , E.D. Earle 14 , 
S.R. Elliott 9 - 19 , H.C. Evans 14 , G.T. Ewan 14 , J. Farine 6 ' 4 , H. Fergani 12 , F. Fleurot 6 , 
R.J. Ford 16 ' 14 , J.A. Formaggio 11 ' 19 , N. Gagnon 19 ' 9 ' 7 ' 12 , J.TM. Goon 10 , K. Graham 14 ' 4 , 
E. Guillian 14 , S. Habib 1 , R.L. Halm 3 , A.L. Hallin 1 , E.D. Hallman 6 , P.J. Harvey 14 , 
R. Hazama 19 ' c , W.J. Heintzelman 13 , J. Heise 2 ' 9 ' 14 - d , R.L. Helmer 18 , A. Hime 9 , C. Howard 1 , 

M.A. Howe 19 , M. Huang 17 - 6 , B. Jamieson 2 , N.A. Jelley 12 , K.J. Keeter 16 , J.R. Klein 17 - 13 , 
L.L. Kormos 14 , M. Kos 14 , C. Kraus 14 , C.B. Krauss 1 , T. Kutter 10 , C.C.M. Kyba 13 , J. Law 5 , 
I.T. Lawson 16 - 5 , K.T. Lesko 7 , J.R. Leslie 14 , I. Levine 4 - 6 , J.C. Loach 12 - 7 , R. MacLellan 14 , 
S. Majerus 12 , H.B. Mak 14 , J. Maneira 8 , R. Martin 14 - 7 , N. McCauley 13 - 12 ^, 
A.B. McDonald 14 , S. McGee 19 , M.L. Miller 11 '*, B. Monreal 11 '' 1 , J. Monroe 11 , 
B. Morissette 16 , B.G. Nickel 5 , A.J. Noble 14 - 4 , H.M. O'Keeffe 12 - 4 , N.S. Oblath 19 , CD. Orebi 

Gann 12 - 13 , S.M. Oser 2 , R.A. Ott 11 , S.J.M. Peeters 12j , A.W.P. Poon 7 , G. Prior 7 < fc , 
S.D. Reitzner 5 , K. Rielage 9 ' 19 , B.C. Robertson 14 , R.G.H. Robertson 19 , M.H. Schwendener 6 , 
J.A. Secrest 13 -', S.R. Seibert 17 - 9 , O. Simard 4 , D. Sinclair 4 - 18 , P. Skensved 14 , T.J. Sonley 11 -™, 
L.C. Stonehill 9 - 19 , G. Tesic 4 , N. Tolich 19 , T. Tsui 2 , CD. Tunnell 17 , R. Van Berg 13 , 
B.A. VanDevender 19 , C.J. Virtue 6 , B.L. Wall 19 , D. Waller 4 , H. Wan Chan Tseung 12 - 19 , 
D.L. Wark 15 - n , P.J.S. Watson 4 , N. West 12 , J.F. Wilkerson 19 - , J.R. Wilson 12 - p , 
J.M. Wouters 9 , A. Wright 14 , M. Yeh 3 , F. Zhang 4 , and K. Zuber 12 - 9 



- 2 - 



1 Department of Physics, University of Alberta, Edmonton, Alberta, T6G 2R3, Canada 
2 Department of Physics and Astronomy, University of British Columbia, Vancouver, BC 
V6T 1Z1, Canada 

3 Chemistry Department, Brookhaven National Laboratory, Upton, NY 11973-5000 
4 Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ot- 
tawa, Ontario K1S 5B6, Canada 

5 Physics Department, University of Guelph, Guelph, Ontario NIG 2W1, Canada 
6 Department of Physics and Astronomy, Laurentian University, Sudbury, Ontario P3E 
2C6, Canada 

7 Institute for Nuclear and Particle Astrophysics and Nuclear Science Division, Lawrence 
Berkeley National Laboratory, Berkeley, CA 94720 

8 Laborat6rio de Instrumentagao e Fisica Experimental de Particulas, Av. Elias Garcia 
14, 1°, 1000-149 Lisboa, Portugal 

9 Los Alamos National Laboratory, Los Alamos, NM 87545 

10 Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 
70803 

11 Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, MA 
02139 

12 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, 
Oxford OX1 3RH, UK 

13 Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 
19104-6396 

14 Department of Physics, Queen's University, Kingston, Ontario K7L 3N6, Canada 
15 Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, UK 
16 SNOLAB, Sudbury, ON P3Y 1M3, Canada 

17 Department of Physics, University of Texas at Austin, Austin, TX 78712-0264 



- 3- 



Received ; accepted 

18 TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada 

19 Center for Experimental Nuclear Physics and Astrophysics, and Department of Physics, 

University of Washington, Seattle, WA 98195 

a Present address: Center for Astrophysics and Space Astronomy, University of Colorado, 

Boulder, CO 

b Present address: Department of Physics, University of Chicago, Chicago, IL 
c Present address: Department of Physics, Hiroshima University, Hiroshima, Japan 
d Present adress: Sanford Laboratory at Homestake, Lead, SD 

e Present address: Department of Physics and Astronomy, Indiana University, South 
Bend, IN 

f Present address: Department of Physics, University of Liverpool, Liverpool, UK 
g Present address: Center for Experimental Nuclear Physics and Astrophysics, and De- 
partment of Physics, University of Washington, Seattle, WA 

h Present address: Dept. of Physics, University of California, Santa Barbara, CA 
"Present address: Department of Physics, Queen's University, Kingston, Ontario, Canada 
jPresent address: Department of Physics and Astronomy, University of Sussex, Brighton, 

UK 

k Present address: CERN, Geneva, Switzerland 

'Present address: Department of Chemistry and Physics, Armstrong Atlantic State Uni- 
versity, Savannah, GA 

m Present address: Dept. of Physics, University of Utah, Salt Lake City, UT 
"Additional address: Imperial College, London, UK 

°Present address: Department of Physics, University of North Carolina, Chapel Hill, NC 

p Present address: Dept. of Physics, Queen Mary University, London, UK 

q Present address: Institut fur Kern- und Teilchenphysik, Technische Universitat Dresden, 



-4- 



SNO Collaboration 



Dresden, Germany 



- 5 - 



ABSTRACT 



We have performed three searches for high-frequency signals in the solar neu- 
trino flux measured by the Sudbury Neutrino Observatory (SNO), motivated by 
the possibility that solar g-mode oscillations could affect the production or prop- 
agation of solar 8 B neutrinos. The first search looked for any significant peak 
in the frequency range 1/day to 144/day, with a sensitivity to sinusoidal signals 
with amplitudes of 12% or greater. The second search focused on regions in 
which g-mode signals have been claimed by experiments aboard the SoHO satel- 
lite, and was sensitive to signals with amplitudes of 10% or greater. The third 
search looked for extra power across the entire frequency band. No statistically 
significant signal was detected in any of the three searches. 



1. Introduction 



Neutrinos are the only w ay known to directly probe the dynamics of the solar 



core (Bahcall and Ulrich 



effect (IMikheev and Smirnov 



1988) , and, through the M ikheev-Smirnov-Wolfenstein (MSW) 
3, 



1986 



Wolfensteii 



19771 ) , they can even carry information about 



the rest of the solar envelope. To date, however, converting measu rements of solar neutri no 



fluxe s into constraints on solar models has proven to be difficult (iBandyopadhyay et al 



20071 ). because of the large number of co- varying parameters upon which such models are 
built. 

A relatively simple signal that could tell us something new about the Sun would be 
a time-variation in the neutrino fluxes. Over the past forty years, measurements made 
by solar neutrino experiments have therefore been the focus of many studies, ranging 
from attempted correlations with the solar sunspot cycle to open searches for signals 



Super-Kamiokande Collaboration 



with periods of weeks or months (ISturrockl 120 03. 2004; Caldwe ll and Sturrockl 12005 



2003 



SNQ Collaboration 



2005al ). The shortest period 



examined to date is roughly one day, where the MSW effect predicts that neutrinos 
propagating through the Earth's core during the night will undergo flavor transformation 
in much the same way they do in the Sun, resulting in a net gain in the flux of electron 
neutrinos (u e s). Although there have been occasional claims of signals on timescales similar 
to known variations in the solar magnetic field, in all cases there have been conflicting 
measurements that show the signals to be spurious or absent entirely. 

We present in this article the results of a search in a new frequency regime for solar 
time variations. Our focus has been on signals whose periods range from 24 hours down 
to 10 minutes. The motivation for such a high frequency search is in part the expectation 
for solar helioseismological variations on scales of order an hour or less, in particular solar 



'gravity modes' (g-modes) ( jChristensen-Dalsgaard 



20031 ). These g- modes are non-radial 



oscillations that are predicted to be confined to the solar core, and thus could in principle 
affect either neutrino production or neutrino propagation. The neutrinos that SNO detects, 
those from 8 B decay within the Sun, are particularly well-suited for our search because they 
are created very deep within the solar core and because their propagation is known to be 
sensitive to variations in the solar density profile through the MSW effect. 



Th e effects of q-modes on so 



Kumar (IBahcall and Kumar 



ar neutrino fluxes have been examined by Bahcall and 



1993), who sought to determine whether g-mode effects could 



explain the apparent solar neutrino deficit, finding that a ny effect was 



Burgess et al. 



ar too small to 



fl2003h looked at 



account for the roughly 60% discrepancy. More recently, 
ways in which a broad spectrum of g-modes could alter the expectation for a solar neutrino 
spectral distortion caused by the MSW effect. Nevertheless, there are at this time no 



explicit predictions as to whether g-modes or any other short-timescale variations could 
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lead to measurable solar neutrino flux variations. 



2. Sudbury Neutrino Observatory 



SNO was an imaging Cherenk ov detector using heavy water (D2O) as both the 



interaction and detection medium (IBoger et al. 



20001 ). The SNO cavern is located in Vale 



Inco's Creighton Mine, at 46°28'30" N latitude, 81°12'04" W longitude. The detector 
resided 1783 m below sea level with an overburden of 5890 ± 94 m water equivalent, deep 
enough that the rate of cosmic-ray muons passing through the entire active volume was just 
three per hour. 

One thousand metric tons of heavy water were contained in a 12 m diameter 
transparent acrylic vessel (AV). Cherenkov light produced by neutrino interactions and 
radioactive backgrounds was detected by an array of 9456 Hamamatsu model R1408 
20-cm photomultiplier tubes (PMTs), supported by a stainless steel geodesic sphere (the 
PMT support structure or PSUP). Each PMT was surrounded by a light concentrator 
(a 'reflector'), which increased the light collection to nearly 55%. Over seven kilotonnes 
(7000 kg) of light water shielded the heavy water from external radioactive backgrounds: 
1.7 ktonne between the acrylic vessel and the PMT support sphere, and 5.7 ktonne between 
the PMT support sphere and the surrounding rock. The 5.7 ktonne of light water outside 
the PMT support sphere were viewed by 91 outward-facing 20-cm PMTs that were used for 
identification of cosmic-ray muons. 

The detector was equipped with a versatile calibration deployment system which could 
place radioactive and optical sources over a large range of the x-z and y-z planes in the 
AV. In addition, periodic 'spikes' of short-lived radioactivity (such as 222 Rn) were added to 
both the light water and heavy water and distributed throughout their volumes to act as 
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distributed calibration sources. 

SNO detected neutrinos through three different processes: 

v x + e~ — ► v x + e~ (ES) 
u e + d^p + p + e~ (CC) 
v x + d^p + n + v x (NC) 



where v x represents u e , or v T . For both the Elastic Scattering (ES) and Charged 
Current (CC) reactions, the recoil electrons were observed directly by their production 
of Cherenkov light. For the Neutral Current (NC) reaction, the neutrons were not seen 
directly, but were detected when they captured on another nucleus. In SNO Phase I (the 
"D 2 phase"), the neutrons captured on the deuterons present within the SNO heavy 
water. The capture on deuterium releases a 6.25 MeV 7 ray, and it is the Cherenkov light 
of the secondary Compton electrons or e + e~ pairs which was detected. In Phase II, (the 
"salt phase"), 2 tonnes of NaCl were added to the heavy water, and the neutrons captured 
predominantly on 35 C1 nuclei. Chlorine has a much larger capture cross section (resulting in 
a higher detection efficiency) for the neutrons. The capture on chlorine also yields multiple 
7s instead of the single 7 from the pure D 2 phase, which aids in the identification of 
neutron events. 

Figure [1] shows the incident 8 B spectrum of neutrinos from the Sun (dotted line), along 
with those that are detected by the CC reaction (dashed) and those that are above the 
effective kinetic energy threshold for the resultant electrons in our detector for Phase II, 
T e ff > 5.5 MeV. T e g is the estimated energy assuming an event consisted of a single electron. 

SNO's depth, very low radioactivity levels, and its ability to perform real-time detection 
made it a unique instrument for observing time variations in solar neutrino fluxes, even at 
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Neutrino Energy (MeV) 

Fig. 1. — Monte Carlo simulation of the 8 B solar neutrino energy spectrum. The dotted 
curve is the incident spectrum of neutrinos from the Sun with an arbitrary normalization, 
the dashed curve shows the spectrum of those detected by the CC reaction before any cut 
on the kinetic energy of the created electron, and the solid curve shows the spectrum of 
neutrinos seen after the application of the kinetic energy threshold for Phase II. 
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the high frequencies we examine here. Above an energy threshold of 5 MeV, the rate of 
events from radioactivity, cosmic ray muons and atmospheric neutrinos, was negligible. The 
rate of solar neutrino events above this threshold was roughly 10/day. 



3. Data Sets 



The event selec tion for the data sets is si milar to that used in our lower-frequency 



periodicity analysis (jSNO Collaboration 



2005a). Events were selected inside a reconstructed 



fiducial volume of R < 550 cm and above an effective kinetic energy of T g > 5 MeV (Phase 
I) or T eff > 5.5 MeV (Phase II), and below T cS > 20 MeV. Additional analysis cuts such as 



fiducial vo 



ume and background rejection fo r these data sets have been described in detail 



elsewhere ( ISNO Collaboration 



2007 



2005bJ). SNO Phase I ran between November 2, 1999 



and May 31, 2001, and with detector deadtimes and periods of high radon removed, we 
recorded a total of 306.58 live days, with 2924 candidate neutrino events. SNO Phase II ran 
from July 26, 2001 to August 28, 2003 for a total of 391.71 live days and 4722 candidate 
neutrino events. Of the 2924 candidate events in the Phase I data, 67% are due to CC 
interactions, 20% to NC interactions, and 9% to ES interactions, with the remaining 4% 
due to backgrounds. Variations in the 8 B neutrino production rate itself will affect all three 
neutrino interactions equally, while variations in the electron neutrino survival probability 
dominantly affect only the CC rate. For the Phase II data the 4722 events consist of 45% 
CC events, 42% NC events, 6% ES events, and 7% backgrounds. 

The time for each event was measured with a global positioning system (GPS) clock 
to a resolution of ~ 100 ns, but truncated to 10 ms accuracy for the analysis. The run 
boundary times were determined from the times of the first and last events in each run with 
a precision of ~ 50 ms. The intervals between runs during which SNO was not recording 
solar neutrino events correspond to run transitions, detector maintenance, calibration 
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activities, and any periods when the detector was off. It is also necessary to account for 
deadtime incurred within a run; for example, deadtime due to spallation cuts that remove 
events occurring within 20 seconds after a muon. This is important for a high-frequency 
periodicity search, as the frequency of occurrence of these deadtimes can approach the scale 
of interest of our search. Therefore, both the run boundaries, as well as the smaller, discrete 
breaks in time due to removal of short-lived backgrounds such as spallation products, define 
the time exposure of the data set, which itself may induce frequency components that could 
affect a periodicity analysis. 



4. Rayleigh Power Approach 



The low-frequency searches for periodicities that have been done by ourselves 
and others typically group the neutrino time series in bins of one to several 
days, and then perform the analysis with methods such as the Lomb-Scargle tech- 



study ( 



nique ( Super-Kamiokande Collaboration 



SNQ Collaboration 



2003 



Sturrock 



2004 ) . In our own low-frequency 



2005al ). we also used an unbinned maximum likelihood technique, 



fitting the time series with periodic functions of varying frequencies, phases, and amplitudes, 
allowing for the detector deadtimes that occurred during calibration runs, power outages, 
and detector maintenance. 

For this high-frequency study, we chose to use an unbinned Rayleigh power approach. 
The Rayleigh power of a time series for a given frequency v is defined as 



z[v) = 



U{u) 



: [(Ej cos I'Kvtif' + (Ej sin 1-nvti 



N N l 

where N is the total number of events in the time series. 

The great advantage of the Rayleigh power approach is its speed, as it requires far 
fewer function evaluations than other unbinned methods like the maximum likelihood 
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technique described above. For this analysis, speed is critical, because to ensure that we do 
not miss a signal we use 1.6 million equally-spaced frequencies spanning a range from 1/day 
to 144/day (one cycle per 10 minutes). To avoid the possibility of a signal falling between 
our sampled frequencies, and thus being missed, the minimum gap between our sampled 
frequencies must correspond to two signals that just decorrelate over the course of SNO's 
running period. With this criterion, the minimum number of frequencies needed for our 
data set in our region of interest is 400,000, and our choice of 1.6 million frequencies was 
made to provide a small degree of oversampling. 

For a time series in which the phase coverage is uniform, the distribution of Rayleigh 
powers for any given frequency follows e~ z , and thus confidence intervals can be easily 
calculated. For the SNO data set, there are significant deadtime intervals, whose durations 
range from months to milliseconds. The sources of these deadtime intervals include the 
period between the Phase I and Phase II data sets (several months), detector calibration 
runs (typically hours to days), power outages (of order one day), maintenance periods 
(hours), and offline veto periods (15 ms to 20 s) imposed to remove events associated 
with the passage of cosmic-ray muons through the detector, interactions of atmospheric 
neutrinos, or bursts of instrumental activity. 

The deadtime structure of the SNO time series means that not all phases of the 
Rayleigh power are equally likely, and thus leads to additional Rayleigh power that is not 
associated with any neutrino signal. Quasi-periodic deadtimes (like those associated with 
calibration and maintenance) can also lead to peaks in the Rayleigh power spectrum. To 
calculate confidence intervals in order to determine the significance of any peaks observed 
in the Rayleigh power spectrum, we must account for these known regions of non-uniform 
phase coverage. 

We have developed an analytic model for the Rayleigh power at a given frequency by 
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treating the Rayleigh power series as a two-dimensional random walk. Each detected event 
is treated as a step in the random walk, with components X = cos 2ixvt and Y = sin 2-nvt. 
For the case of uniform phase coverage, the central limit theorem implies that for a large 
number of steps (N) the distribution of final positions will be given by a two-dimensional 
Gaussian, whose means, variances, and covariance are: 



^ x = 2tt ' ^ COS( t ) = Q ( 2 ) 



1 f 2w . , , . 1 



and 

a * = 2^ J d ^ cos2 ^ = i ^ 

1 f 2w 

cov(x, y) = — / d(p cos sin = (4) 
2tt J 

leading to a simple distribution of final positions given by 

f(X,Y) = — i exp(-X 2 /2iV^) exp(-r 2 /2iV^) 

Z7TiV O x O y 
= 1 C -(X 2 +Y 2 )/N 



Nti 

1 

Ntt 



e -U 2 /N 



The distribution of z = U 2 /N can then be obtained by integrating over all values of X 
and Y which satisfy z < U 2 /N < z + dz, by changing variables from X and Y to ip and U : 

f(X, Y) dX dY -> /([/, V) rfC/ ^ = 4re- u2,N U dU dip. (5) 

N71 

Integrating this expression over dip from to 2ti and changing variables from U to z = U 2 /N 
gives f(z) dz = e~ z dz, which is the simple exponential distribution expected for the 
Rayleigh power distribution. 

For the case of non-uniform coverage, the means and covariance of the distribution are 
no longer simple. If we call the normalized phase-weighting function g((f>), where g((f>) = 1 
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for uniform phase coverage, then we have 



fJ-x = — d<j) g{4>) cos0 (6) 
2 ?r Jo 



1 f 2w 

H y = — J d<p #(0)sin0 (7) 



^ r 2it 



^ 2 = — y d(j> g{<t>) (cos (f)-fi x ) 2 (8) 



Oy = — j d(f) g ((f)) (sin - /i y ) 2 (9) 



1 r 71 

cov(x,y) = —J d(j) £f(0)(cos0 - /^)(sin0 - /i y ) (10) 

The function g(0) is determined by the detector's deadtime window, with <p = out. The 
mean for X, for example, is given by 



runs „i 



stop,j 



H x = - / dtcosut (11) 

j=l ^ ^ start, j 



where T is the total livetime, and the sum is over all data taking runs in the data set, 
integrating from the start to the stop time of each run. 

The Rayleigh power distribution for the non-uniform phase coverage case is proportional 
to e~ x / 2 , and x 2 f° r the Rayleigh power distribution can be written as 

( X — N/jl x , 

X 2 (X, Y) — (X — Nn x , Y - Nfi y ) V xy - X | (12) 

\Y-Nfi y 

The inverse of the covariance matrix is given by 
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v - 1 



-1 

Na x 2 Ncov(x, y) 



Ncov(x, y) Na, 



2 



making our x 



2 



2 _ (X - Nfx x ) 2 a y 2 + (Y- Nfx y )W - 2(X - Nfx x )(Y - iV^)cov(x, y) 
X Ncr x 2 a y 2 -Ncov 2 (x,y) 1 ' 

Transforming into our integration variables z = U 2 /N and i[) gives the rather unwieldy 
probability density function for the Rayleigh power at a given frequency: 

f(z)dz = ^dz f\H<*m/m*A<xm/m^m-{<*m/m* dlf) (14) 

C Jq 

where 

a.\(ij)) = a y 2 cos 2 ip + a x 2 sm 2 ip — 2 cosip smipcov(x,y) 
&2{ip) = — 2(<Tj / 2 /i a; cos ij) + <J x 2 \i y sin^ + [ji x sin^ + n y cos^)cov(x, y)) 
a 3 (^) = cr y 2 [i x 2 + a 2 [iy 2 + /i x /iyCov(x, y) 

and C is a normalization constant. In this formalism the effects of detector deadtime are 
entirely accounted for through the means and covariance matrix for the variables X and Y. 
In the analysis of data from the two combined SNO data-taking phases, the difference in 
event rates between Phase I and Phase II is accounted for by separating the terms of the 
analytic form according to phase, or: 

Nfl x ,y — > N D2 o(fi x ,y) D2 q + Nsalt(Vx,y) Sa i t 



(and similarly for all variance and covariance terms). Here, both g(<p) and N have been 
separated according to phase, effectivel y introduc i ng a r ate-dependent weighting factor. For 
further details on this method, see Ref. lAnthonyl (120081 ). 
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Figure [2] shows the function f(z) describing the Rayleigh power distributions for two 
different frequencies. The distributions were generated with a Monte Carlo simulation 
including the SNO detector's full deadtime window. For the plot on the left, the deadtime 
contributes a significant amount of power due to the periodicity in SNO's operations 
schedule, while for the high frequency bin on the right the deadtime does not change the 
function much from its simple e~ z distribution. 

To determine a specific confidence level, CL, for a given observed Rayleigh power, z , 
we solve the equation 

CL= / f(z)dz (15) 
Jo 

As a test of our analytic model, we have calculated the confidence levels for all 1.6 
million frequencies of a Monte Carlo simulation that includes the full SNO deadtime 
window. Figure [3] shows the distribution of these confidence levels, which is gratifyingly 
flat with a mean that is 0.50, thus showing that the analytic random walk model correctly 
distributes the confidence levels across the whole Rayleigh power spectrum. 



5. Open Single Peak Search 

Our first search looked for a significant peak at any frequency in our Rayleigh power 
spectrum. While Equation [15] gives the confidence level at any specific frequency, when 
testing our 1.6 million sampled frequencies there is a substantial trials penalty, making it 
exceedingly likely that at least one of the frequencies will by chance have an apparently 
large power. To determine this penalty exactly, we would need to know how many of our 
1.6 million sampled frequencies are independent, which is a complex task. 

Instead, to address this trials penalty we use 10,000 null-hypothesis Monte Carlo 
simulations to determine the probability of observing a statistically significant peak at 
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any of the 1.6 million sampled frequency in the absence of a true signal. For a given 
null-hypothesis Monte Carlo simulation, we assign a frequency-specific confidence level to 
each sampled frequency according to the prescription given above in Section H] (see Eq. [T5|) . 
Then for each null-hypothesis Monte Carlo simulation, we record the peak that has the 
highest confidence level, and then plot the resultant distribution of these highest-peak 
confidence levels for all 10,000 null-hypothesis Monte Carlo simulations. The resultant 
distribution of confidence levels is shown in Figure HI As seen from this figure, virtually 
every simulation yields at least one peak with an apparent significance of at least 99.999%, 
just by chance. To determine the true data-wide confidence level, taking into account our 
entire sample of 1.6 million frequencies, we place a cut on the distribution of Figure H] 
that corresponds to our desired trials-weighted (data-wide, rather than frequency-specific) 
confidence level for a significant signal. In Figure HI the cut shown corresponds to the 
frequency-specific confidence level needed by the maximum peak in a power spectrum in 
order for it to be above the 90% CL detection threshold. In other words, for a data set that 
contains no periodicity at any frequency, there is a less than 10% probability that the most 
significant individual peak will have a frequency-specific confidence level in excess of this 
cut value. 

Figures [5] through [10] show the Rayleigh power spectrum for the combined SNO Phase 
I and Phase II data sets, broken up into six segments each corresponding to roughly 267,000 
sampled frequencies. 

Figure [TT] shows the peak with the highest confidence level, and the corresponding 
threshold for that peak to be above the data- wide 90% CL to be considered significant. The 
data- wide CL of this peak is only 2%, thus we see no evidence of a significant peak in our 
data set. 

To determine our sensitivity to a signal, we ran Monte Carlo simulations with fake 
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sinusoidal signals, of form 1 + A sin out, of increasing amplitude, looking for the point at 
which our method would claim a discovery. In Figure [Tj2] we show our sensitivity for two 
criteria: the amplitude required to make a 99%-CL discovery 90% of the time, and the 
amplitude for a 99%-CL discovery 50% of the time. We are substantially limited in this 
open search by the trials penalty; we need a signal of 12% amplitude to make a 99% CL 
detection 90% of the time. The bands shown in Figure [T21 indicate the degree of variation 
among frequencies of the sensitivity, which is affected by the underlying power spectrum in 
each bin as discussed above in Section HI 



Directed Peak Search 



There have been recent claims by the GOLF/SoHO collaboration of possible 



signa tures of q-mode oscillations, based on 



sets ( jGarcia et al 



2001 



Gabriel et al. 



2002 



analyses of long-term he 



Turck-Chieze et al. 



2004 



i oseismological data 



Mathur et al 



2003), 



20081 ). 



as well as supporting claims by the VIRGO/SoHO collaboration tlGarcia et al. 
Looking for such specific signals in our data set using our Rayleigh power approach has an 
advantage in that we no longer need the 1.6 million frequencies used above, but rather can 
look in a narrow window that has a smaller trials p enalty. We have thus take n a narrow 



band around the reported persistent GOLF signals (I Jimenez and Garcia! 120081 ). from 18.5 
to 19.5/day (roughly 214 to 225 /iHz), and have repeated our analysis. We again find no 
significant signal, with the highest peak having a trials-weighted CL of just 58%. The 
Rayleigh power spectrum around this peak, as well as its superimposed 90% CL are shown 
in Figure [131 Figure dH shows our sensitivity plots for this directed search, which are 
slightly better than in Section [5] because of the reduced trials penalty. We conclude that if 
the detection claimed by SoHO is in fact evidence of a g-mode, the effect of this particular 



mode of oscillation on the neutrino flux is less than 10% amplitude variation, at 99% CL. 
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7. Broadband Search 

The two searches described above require that the signal be predominantly sinusoidal 
and monotonic. It is possible that high freq uency behavio r in th e Sun spans a large band 



of frequencies, and in fact may be 'noisy' 



Burgess et al 



(120031 ) have investigated how 



such noise might affect the neutrino survival probabilities within the Sun due to the matter 
or MSW effect. We have therefore looked at the distribution of confidence levels across 
our entire range of 1.6 million frequencies. Like Figure [31 we expect that in the case of 
no broadband signal the distribution of confidence levels will be flat, with a mean of 0.50. 
Figure [15] shows this distribution now for our combined SNO Phase I and Phase II data 
sets. As can be seen clearly in the figure, the distribution is flat, with a mean very close to 
the 0.5 expected. As a comparison case, we show in Figure [TBI what the confidence level 
distribution of a 'noisy' sun would look like, for several different amplitudes of Gaussian 
white noise. Our white noise model is shown in Figure [TT| for the lowest amplitude (0.1%) 
shown in Figure [161 We have spread the noise across 400,000 independent frequencies, 
roughly the number of independent frequencies we expect in the power spectrum. The 
rms noise power from our model is 0.45, in units of SNO's measured total 8 B neutrino 
flux (~ 5 x 10 6 z/ cm~~ 2 s~ 1 ). In terms of power per unit bandwidth, this corresponds to 
~ 6 x 10 7 u cm~ 2 s~ 1 Hz~ 1//2 . As is evident in Figure [151 t ne distribution of confidence levels 
in the data is consistent with no distortion of signal due to noise. 



8. Conclusions 

We have performed three searches for high-frequency signals in the 8 B solar neutrino 
flux, applying a Rayleigh Power technique to data from the first two phases of the Sudbury 
Neutrino Observatory. Our first search looked for any significant peak in a Rayleigh Power 
spectrum from frequencies ranging from 1/day to 144/day. To account for SNO's deadtime 
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window, we calculated the expected distribution of power in each bin of the Rayleigh Power 
spectrum using a random walk model, thus allowing us to assign confidence levels to the 
observed powers. We found no significant peaks in the data set. For this 'open' peak search, 
we had a 90% probability of making a 99% CL detection of a signal with an amplitude of 
12% or greater, relative to SNO's time-averaged neutrino flux. 

In a second search, we narrowed our frequency band to focus on a region in which 
g-mode signals have been claimed by experiments aboard the SoHO satellite. The examined 
frequency range extended from 18.5/day to 19.5/day. Again, no significant peaks in the 
Rayleigh Power spectrum were found, and our sensitivity for this 'directed' search gave us 
a 90% probability of making a 99% CL detection for signals whose amplitudes were 10% or 
larger, relative to SNO's time-averaged neutrino flux. 

Our third search examined the entire range of frequencies from 1/day to 144/day, 
looking for any evidence that additional power was present across the entire high-frequency 
band. To do this, we used the distribution of frequency-specific confidence levels, determined 
using our random walk model. We found that, as expected for no high-frequency variations, 
this distribution was flat. We showed that for a simple Gaussian white noise model, the 
confidence level distribution would be noticeably distorted even when the amplitudes of the 
contributing frequencies have an rms as small as 0.1%. 

This research was supported by: Canada: Natural Sciences and Engineering Research 
Council, Industry Canada, National Research Council, Northern Ontario Heritage Fund, 
Atomic Energy of Canada, Ltd., Ontario Power Generation, High Performance Computing 
Virtual Laboratory, Canada Foundation for Innovation; US: Dept. of Energy, National 
Energy Research Scientific Computing Center; UK: Science and Technologies Facilities 
Council. We thank the SNO technical staff for their strong contributions. We thank Vale 
Inco, Ltd. for hosting this project. 
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(a) Distribution of powers at sampled fre- (b) Distribution of powers at sampled fre- 

quency = 1.000089 day" 1 . quency = 1.046359 day" 1 . 

Fig. 2. — Distribution of specific frequencies' Rayleigh powers in SNO Monte Carlo data sets 
(black) vs predictive analytic form (red line). 
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Fig. 3. — Distribution of Rayleigh power confidence levels for all frequencies of a Monte 
Carlo simulation of the SNO data set. 
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Fig. 4. — Distribution of maximum confidence levels from Rayleigh analysis of 10,000 signal- 
free Monte Carlo data sets. By building this distribution of maximum confidence levels, 
we can determine a 'confidence level of confidence levels' and account for the trials penalty 
in our generation of the data-wide confidence level. The frequency-specific confidence level 
corresponding to the data-wide confidence level of 90% is shown with the superimposed 
dashed red line. 
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Fig. 5. — First of six figures. Rayleigh power spectra for first subsection (frequencies between 
1 and 15 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 day" 1 ), for 
combined SNO Phase I and Phase II data sets. The entire range has been broken down into 
six individual frames for easier inspection, with the first frame slightly more zoomed-in due 
to the presence of more underlying activity in this region of lower frequencies. The black 
line indicates data, and the upper red line designates the level at which a detection would 
have a confidence level of 90%. The peaks at low frequencies, specifically those at lday -1 
and 2day _1 , represent SNO-specific periodicities due to daily run-taking schedules. (This 
is clearly not evidence of a signal, as the red CL=90% line, generated from null-hypothesis 
Monte Carlo, also follows these peaks). 
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Fig. 6. — Second of six figures. Rayleigh power spectra for second subsection (frequencies 
between 15 and 41 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 day -1 ), 
for combined SNO Phase I and Phase II data sets. The black line indicates data, and the 
upper red line designates the level at which a detection would have a confidence level of 90%. 
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Fig. 7. — Third of six figures. Rayleigh power spectra for third subsection (frequencies 
between 41 and 67 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 
day -1 ), for combined SNO Phase I and Phase II data sets. The black line indicates data, 
and the upper red line designates the level at which a detection would have a confidence 
level of 90%. 
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Fig. 8. — Fourth of six figures. Rayleigh power spectra for fourth subsection (frequencies 
between 67 and 93 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 day -1 ), 
for combined SNO Phase I and Phase II data sets. The black line indicates data, and the 
upper red line designates the level at which a detection would have a confidence level of 90%. 
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Fig. 9. — Fifth of six figures. Rayleigh power spectra for fifth subsection (frequencies between 
93 and 119 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 day -1 ), for 
combined SNO Phase I and Phase II data sets. The black line indicates data, and the upper 
red line designates the level at which a detection would have a confidence level of 90%. 



-28- 



Si 40 




Frequency (day 1 ) 



Fig. 10. — Last of six figures. Rayleigh power spectra for sixth subsection (frequencies 
between 119 and 144 day -1 ) of entire range of frequencies sampled (from 1 day -1 to 144 
day -1 ), for combined SNO Phase I and Phase II data sets. The black line indicates data, 
and the upper red line designates the level at which a detection would have a confidence 
level of 90%. 
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Fig. 11. — Zoomed- in region of the Rayleigh power spectrum for the highest-significance peak 
in the SNO data set, which was detected at frequency=103.384 day' 1 , with a confidence level 
of 2%. The horizontal red line indicates the frequency-specific powers needed for a peak to 
be above the data-wide 90% CL. 
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Fig. 12. — SNO's sensitivities to a high-frequency periodic signal in the combined data sets 
(Phase I, or D 2 0, and Phase II, or salt) for the entire frequency search region. The cyan 
band shows the calculated sensitivity at which we detect a signal 50% of the time, with 99% 
CL, and the red band shows the calculated sensitivity at which we detect a signal 90% of 
the time, with 99% CL. The width of the bands represents the range of variation of the 
sensitivity, which varies rapidly with frequency, across the frequency regime. 
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Fig. 13. — Rayleigh power spectrum for 'directed' high-frequency search, in black. The line 
corresponding to detection with 90% CL is shown in red. The highest peak in the power 
spectrum is found at a frequency of 19.2579/day with a CL of 58%. 
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Fig. 14. — SNO's sensitivities to a high-frequency periodic signal in the combined data sets 
(Phase I, or D2O, and Phase II, or salt) for the directed high frequency search region. The 
cyan band shows the calculated sensitivity at which we detect a signal 50% of the time, with 
99% CL, and the red band shows the calculated sensitivity at which we detect a signal 90% 
of the time, with 99% CL. The width of the bands represents the range of variation of the 
sensitivity, which varies rapidly with frequency, across the frequency regime. 
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Fig. 15. — Distribution of confidence levels for all 1.6 million frequencies, for the SNO 
combined-phase (D 2 and salt) data set. 
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Fig. 16. — Distribution of confidence levels for all 1.6 million frequencies in a SNO white- 
noise Monte Carlo, with several signal amplitudes. The SNO combined-phase (D2O and 
salt) data confidence level distribution is shown in black for comparison. 
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Fig. 17. — Time domain plot of our Gaussian white noise model, with amplitude of 0.1%, 
distributed over 400,000 frequencies. Only ten days are shown here. The rms noise power 
for the model shown here is 0.45 in units of SNO's measured total 8 B neutrino flux. 
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